Measurement of Counting Statistics of Electron Transport in a Tunnel Junction 
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We present measurements of the time-dependent fluctuations in electrical current in a voltage- 
biased tunnel junction. We were able to simultaneously extract the first three moments of the tunnel 
current counting statistics. Detailed comparison of the second and the third moment reveals that 
counting statistics is accurately described by the Poissonian distribution expected for spontaneous 
current fluctuations due to electron charge discreteness, realized in tunneling transport at negligible 
coupling to environment. 



Photon counting statistics Q, 0] is a key technique of 
quantum optics, which is routinely used to study tempo- 
ral and spectral distribution of electromagnetic field, and 
to characterize the complexity of optical states, such as 
photon coherence, entanglement, and squeezing. In con- 
trast, the electron counting and noise statistics, which are 
expected to provide new insights into quantum transport, 
is essentially at infancy as experimental subject, with the 
first advances made only recently 0, 0, IE S ■ 

Electron counting proves to be far more challeng- 
ing than that of photons primarily because of the ex- 
tremely high frequency of electron passage events at typi- 
cal current intensity, which requires fast charge detection. 
While in some cases Coulomb blockade in quantum dots 
can be used to localize electrons and suppress tunnel- 
ing, bringing the tunneling frequency down to the radio 
frequency range 0, 0, , the time resolution needed to 
measure free, nonlocalized electrons remains a challenge. 

Another difficulty stems from the simple fact that, 
while photons do not interact, the electrons do. The elec- 
tric field fluctuations produced by the electrons which are 
being measured can propagate out to other parts of the 
circuit ('environment') and perturb it. In turn, the noise 
due to the environment, modulated by the signal, can 
couple back to the region of interest, strongly affecting 
the measured signal % H IS GH 111 G3 

The electron problem is especially rich and intriguing: 
unlike photons, the current fluctuations are measured 
without extracting electrons out of the system. The elec- 
trons remain part of the many-body system while being 
detected, allowing quantum phenomena to fully manifest 
themselves in electric noise. This leads to a number of 
dramatic effects observed in electric noise, such as, no- 
tably, elementary charge transmutation in the Fractional 
Quantum Hall effect [l]3, HH 113 an d charge doubling in 
the Andreev scattering regime in NS junctions |l6|. 

The regime in which electric fluctuations occur spon- 
taneously due to charge discreteness, rather than due to 
thermal fluctuations, is realized at sufficiently low tem- 
peratures [13, llo, LLi| • It was first demonstrated about 
10 years ago in semiconductor point contacts [13, El 
and in mesoscopic wires |22l Current fluctua- 



tions in this regime are traditionally analyzed using the 
noise frequenc y sp ectrum. However, a more detailed 
description[3 l25| is provided by the statistics of the 
charge q(r) transmitted through a conductor during a 
fixed time interval r which, in principle, can be small 
or large compared to the time tq = e/I. The low fre- 
quency shot noise is just the second central moment of 
the transmitted charge distribution P(q) at long times 
T > t : 



S {2) = (1*1^=0 



<Ag 2 (r)) _ ((q(r)-qr) 



(1) 



with (...) a shorthand for J ...P(q,r)dq. The counting 
distribution P(q), which encodes the full information 
about noise statistics, was studied in various regimes, in- 
cluding mesoscopic scattering |24j , photon-assisted trans- 
port [2g, and transport in NS junctions [Hl^g, where 
charge doubling takes place. Specifics of the tunneling 
regime were considered in |29|, |3fjj| . 

The present work extends the shot noise measurements 
beyond the second moment in a way uncorrupted by 
the presence of electromagnetic environment. The re- 
sults agree with the expectation for Poisson process, thus 
paving the way to investigation of counting statistics. We 
measure current fluctuations in a tunnel junction by de- 
tecting the probability distribution of voltage across a 
load resistor, P(V). The latter is directly related to P(q) 
when the load resistance is made much smaller than the 
tunnel resistance. The knowledge of P(q) allows, in prin- 
ciple, to obtain all moments of q. However, with the mea- 
surement times t ^> ro, due to the central limit theorem, 
the high moments become increasingly dominated by the 
lower order moments. This makes the irreducible parts of 
the moments, the cumulants, which contain new informa- 
tion, increasingly difficult to extract. Thus here we shall 
focus on the third central moment (Ag 3 ) = ((q — q) 3 ) 
which coincides with the third cumulant ((q 3 ))- 

The expected statistics is different in the voltage- 
biased and current-biased regimes. The former case is 
described by the rate of the attempts to pass equal to 
eV/h and a binomial distribution of successes [24j. At 
weak transmission, realized in tunnel junction, the dis- 
tribution assumes Poisson form, with the low frequency 
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spectral density of the fc-th cumulant, corresponding 
to the 'long' measurement time r > To, expressed as 



S (k) 



«<?* 



(2) 



In contrast, in a current-biased sample, the rate of suc- 
cesses is fixed at I/e, and the attempts to pass, described 
by fluctuating time-dependent voltage across the sample, 
are characterized by Pascal distribution A general 
cascade approach to high-order noise statistics in the dif- 
fusive regime was developed in Ref. lal , the role of exter- 
nal circuit was considered in Ref. |lOl | (see also ^H). 

In the first measurement of the third moment S^ 3 ' 
Reulet et al. 0, Q used a low impedance tunnel junc- 
tion (i? s ~ 50 Ohm) as a noise source, with a parallel 
50 Ohm load impedance Ri of the cable used to feed the 
voltage to an external amplifier. The initial results Q 
distinct difference from the Poissonian, in both magni- 
tude and sign, suggested that S 1 - 3 ' is dominated by the 
effects extraneous to the junction. A theoretical inves- 
tigation 0, followed which clarified the importance 
of the external circuit for correct interpretation of the 
results, obtaining the third moment of the form 



SO) = - 3R Syl (S^ + S t 



(2 



h B(S. 
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dV 



(3) 



with R St i = R s Ri/(R s + Ri) the impedance of the sam- 

(2) (2) 

pie in parallel with the load, Ss and S} the sample 

(3) 

and load noise, and Ss generated by the voltage-biased 
sample. The load resistor, being macroscopic, is not ex- 
pected to produce a third cumulant [jUl^. The measure- 
ment 0, , which due to Ri ~ R s was neither fully in 
the voltage-, nor in the current-biased regime, was dom- 

(2) 

inated by the voltage-dependent S s in the last term of 
Eq.(|3J|. Only at room temperature, due to low dS^ /dV, 
the result had the sign of Ss 3 ^ . 

To measure the intrinsic S^ free of the admixture 
of the second moment, we use a new method suggested 
and analyzed in Ref.[2!| (Fig|TJ. Current fluctuations 
generated by the sample (voltage-biased tunnel junction 
of high impedance) produce voltage fluctuations on the 
load resistor Ri which are amplified and analyzed with 
computer. The statistics of voltage fluctuations on Ri is 
identical to that of the intrinsic current fluctuations in 
the junction, provided that Ri is much smaller than the 
junction differential resistance. 

The main source of errors in the measurement of the 
fc-th cumulant ({q k }} of the distribution P(q) is statis- 
tical. In order to estimate the signal-to-noise ratio the 
measured value ((q k )} should be compared to its variance 
var(Aq fc ) due to both sample and amplifier noise. The 
variance is expressed through the central moments of the 
order 2k. The variance of an even order for a generic dis- 
tribution can be estimated with the help of the central 
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FIG. 1: Schematic of the experimental setup. The current 
source drives constant average current I while the capacitor 
C v fixes the voltage across the sample and the load at relevant 
frequencies. The bandwidth is determined by the load resis- 
tance Ri and stray capacitance C 3 . The voltage fluctuations 
across the load resistor are amplified and digitized by a 12 bit 
analog-to-digital convertor (A/D) obtained from Ultraview 
Corp. with 5 ns sampling time and 20 ns interval between the 
samples. The amplifier input impedance is determined by Ri . 



limit theorem, using Gaussian statistics: 



var 



{Aq k ) = ({Aq 2k )) 1/2 ~ {(2k - l)!!) i/2 {Aq 2 ) k / 2 (4) 

Here the variance (Ag 2 ) in Eq.(£j} is due to the sample, 
amplifier and load noise, — Sf^ + sj 2) + S^ 2 \ as- 
suming all three to be uncorrelated. The signal-to-noise 
ratio for a single measurement of the third cumulant is 
thus estimated as a ratio of S3T to 15 1 / 2 (S , 2t) 3 / 2 , and 
therefore it is beneficial to reduce the sampling time r to 
improve sensitivity, in accord with the central limit theo- 
rem. However, the amplified signal is correlated at short 
times due to finite bandwidth of the input circuit. This 
makes the effective sampling time restricted from below 
by the sample parasitic capacitance, C s , of both intrin- 
sic and stray kind, and by the effective output resistance 
of the sample in parallel with the load resistance of the 
amplifier, i?/, so that r c g ~ RiC s . Repeating the mea- 
surement N times would further reduce fluctuations by a 
factor y/N, assuming statistically independent successive 
measurements. However, since the measurements sepa- 
rated in time by less than r e g are correlated, the maximal 
improvement is of the order of y/T /t c h, where T is the 
measurement time. In the case of a high impedance tun- 
nel junction, the noise is dominated by the resistor Ri 
thermal noise, 2ksT/Ri. The measurements were per- 
formed at 4.2 K to reduce both this noise and the noise 
of the amplifier. Ignoring for now both the shot noise 
produced by the sample and the amplifier noise, we esti- 
mate the signal-to-noise ratio as 



S/N 



(5) 



/ 15(2fc s T/i? ; ) 3 / 2 V^i 

Replacing r by r e g and plugging in Eq.© we find the 
signal-to-noise ratio S/N oc J?j C" 1 . It is therefore 
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FIG. 2: Measured value of the noise, as a function of 

the current. Different markers represent the data of three 
different measurements. The straight line S 1 ' 2 ' = el, corre- 
sponding to the expected noise intensity |3l|. is shown for 
guidance. Upper inset: raw voltage histogram sample; Lower 
inset: cleaned histogram, illustrating the result of normalizing 
with the A/D converter calibration (see text). 

clear that it pays to decrease C s and increase Ri to im- 
prove S/N. We placed a cold amplifier in the vicinity 
of the sample to reduce as much as possible the capaci- 
tance C s . The choice of the Ri is restricted by the desire 
to keep most of the signal at frequencies well above 1// 
corner of the amplifier. We chose Ri = 9.1kOhm which 
together with the stray capacitance of about 4 pF gives 
r ff = RiC s ~ 40ns and signal bandwidth (27rT e ff) _1 « 
4 MHz. To keep the sample under voltage bias, we in- 
troduced a capacitor C v which kept the voltage constant 
across the sample-load circuit at all relevant frequencies. 

Wc used tunneling through a 30A thick SiC>2 gate oxide 
of a p-channel Si field-effect transistor to produce a shot 
noise. In this system tunneling occurs only under nega- 
tive gate voltage required to induce a hole channel. The 
differential resistance of the barrier R s = (dl / dV)^ 1 > 
10 7 Ohm was much bigger then Ri, placing the sample se- 
curely in the voltage-biased regime. Indeed, the maximal 
contribution of the second term in Eq.(|3J|, 

is estimated as 5.2 x 10 -48 A 3 /Hz 2 , which is two orders 
of magnitude smaller then e 2 I (Fig.[2J). 

We record the probability distribution function P(V a ) 
of the amplified voltage V a (Fig.|2 upper inset). To 
clean the obtained histogram we correct for nonequal bin 
widths of the A/D converter by calibrating the latter 
against linearly swept signal (Fig.|3 lower inset). The 
effect of normalizing P(V a ) with the A/D converter cali- 
bration, illustrated in Fig. flower inset, is two-fold. First, 
the histogram loses the erratic features ('grass') on the 




FIG. 3: Measured third cumulant ' of the transmitted 
charge, obtained separately for different current directions. 
(Markers are the same as in Fig.|21 the straight line is S 1 ' 3 ' = 
e 2 I.) Upper inset: S^ 3 ' vs. I without amplifier nonlinearity 
correction; Lower inset: normalized histogram of the linearly 
swept signal, used to calibrate the A/D converter (see text). 

smallest scale. Simultaneously, the envelope is somewhat 
corrected on large scale. The latter effect, due to aver- 
aging, has longer stability time, which is fortunate, since 
the 'grass' features, while less stable, were found to have 
little effect on the second and the third cumulant. 

For a linear circuit, the amplified voltage V a {uj) and the 
input current i{u>) Fourier components are proportional: 
V a {uj) — Z{oj)K(u))i(uj), with Z the load impedance and 
K the amplification. Therefore the third cumulant is 
related to the respective spectral densities of V a as 

(V®) = j^JJ Ji{w)A{u>')A{-u-w')<L>dM>' (7) 

with A(ui) — Z(lu)K(lu). For the broadband amplifiers 
used, the amplification K(lo) was almost frequency inde- 
pendent between the low ~ 0.5 MHz and high ~ 25 MHz 
frequency cutoffs intentionally introduced to filter the 
1// amplifier noise and the wide-band noise at fre- 
quencies well above (27rT e fi) _1 . The complex product 
Z(uj)K(u>) was obtained from calibration, introducing a 
known current through a small capacitor C c typically of 
2.4 pF. We then used Eq.Q for vi 3 "* and a similar expres- 
sion for K (2) to extract and . The second cumu- 
lant combines the shot, thermal and amplifier noise 
contributions. Since only the former depends on current, 
the amplifier noise and the thermal noise, obtained from 
at zero current, can be easily subtracted [12 The 
resulting noise SW, shown in Fig.0 varies linearly with 
current as expected (eV ~ 4V » fc^T at all currents). 
The measured value agrees with the expected value, el, 
within calibration error of 5%. 

In order to extract properly the third cumulant of the 
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voltage, one should account for the effect of amplification 
nonlinearity which can mix 5 (2) with S {3 \ Let us con- 
sider the amplifier nonlinearity, V out = K(V in + V 2 n /U), 
with 1/U the nonlinearity parameter, which which con- 
tributes to the cumulants (SV 2 ut ) and (SV 3 ut ) as follows: 



(SV 2 ut )/K 2 - (5V£) = 2(5Vl)/U + 3(5V£) J /U J (8) 



(9) 



The right hand side of the Eq. @ would affect the value 
calculated from (5V„ ut ). The parameter U can be 
estimated by applying an ac signal and measuring its 
second harmonic. We took special care to reduce nonlin- 
earity of the amplifier track, especially the last amplifier 
in the chain. We managed to make it comparable to the 
nonlinearity of the cryogenic amplifier which was differ- 
ent for different measurements, depending on the regime. 
We found U to exceed 20 mV, and thus having negligible 
effect on S^ 2 '. However, estimates show that the nonlin- 
earity correction to S^ 3 ' can be as large as 10 -45 A 3 /Hz 2 . 
We attribute the third cumulant at zero current (see the 
upper inset of Fig. 01 to this nonlinearity and determine 
U as the ratio 9(5V 2 n ) 2 I=0 / (SV 3 ut ) I=0 . Using this U we 
subtract the nonlinearity contribution, Eq.JHJ), from the 
data at all currents. This procedure, as illustrated in 
Fig. |3 upper inset, restores zero value at I = 0, but has 
little effect on the slope of the dependence vs. I. We 
therefore conclude that the nonlinearity, while necessary 
to account for to improve accuracy, is not strong enough 
to compromise the measurement of dS^/dl. The re- 
sults, shown in Fig. [31 are found to be in excellent agree- 
ment with the Poissonian third cumulant e 2 I. 

As an additional check, we also reversed the current 
direction. Since we had to apply negative voltage on the 
gate to induce holes, it required rebonding the sample 
(see Fig.|3 upper inset). We found the third cumulant 
to be an odd function of current, as expected. 

Note added: Finally, we note that the literature is not 
entirely unanimous regarding the Poissonian character of 
tunneling transport. Dissent is exemplified by Ref. [jjj] . 
predicting a new "quantum regime" at sufficiently low 
frequencies, limited by inverse flight time in a noninter- 
acting fermion model. Ref. |3^| obtains S^ 3 ' for the tun- 
neling current of sign opposite to Poissonian and 
also of much smaller magnitude, quadratic in transmis- 
sion rather than linear as in Eq.J2J above. The condi- 
tions stated in Ref. |33l | are fulfilled in our experiment: 
the time interval between electron tunneling and its de- 
tection, estimated from EM signal propagation speed, is 
of order 3 • f 0~ n s, which is much shorter than the sam- 
pling time, 5 ns, placing the measurement securel y at low 
frequency in the sense of Ref. (3j| . The results |33| are 
thus not in accordance with our observations. 



In summary, we present the first measurements of the 
charge counting statistics in voltage-biased tunnel junc- 
tion up to the third cumulant. The results, obtained 
by analyzing the distribution of transmitted charge, are 
in excellent agreement with expectations for Poissonian 
process, making electron counting statistics amenable to 
experimental investigation. 
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Science Foundation and (D.S. and Yu. B.) by the Lady 
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